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In-Sltu  Material  Characterization 
For  Pavement  Evaluation 
bv  the 

Spectral -Analysis -of -Surf ace -Waves  (SASW) 

Method 

I .  INTRODUCTION 

A .  Background 

In  situ  material  characterization  is  necessary  for  the  structural 
evaluation  of  existing  pavement  systems.  In  the  past,  the  accepted 
method  of  evaluating  pavements  has  been  through  destructive  testing, 
where  strength  parameters  are  determined  from  field  and  laboratory 
evaluation  of  pavement  layers .  More  recently  a  growing  acceptance  of 
nondestructive  testing  (NDT)  as  a  means  of  characterizing  pavement 
materials  has  led  to  the  development  and  widespread  use  of  various  NDT 
devices  and  pavement  evaluation  procedures.  NDT  offers  the  advantages 
of  causing  less  interference  to  normal  traffic  operations  and  being 
less  costly  than  destructive  testing  procedures. 

In  the  late  1970' s,  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  developed  a  methodology  for  evaluation  of  light  aircraft 
pavements  based  upon  multilayered  elastic  models  and  limiting 
stress/strain  criteria  (Bush,  1980).  Modulus  values  can  be  determined 
using  surface  deflections  obtained  from  a  number  of  commercially 
available  test  devices,  one  of  the  most  popular  being  the  falling 
weight  deflectometer  (FWD) .  This  procedure  is  currently  being  used 
extensively  by  the  military  for  airfield  pavement  evaluation  and  is 
also  being  widely  used  by  many  Universities  and  engineering  firms  for 
pavement  research  and  evaluation.  However,  with  the  layered  elastic 
evaluation  procedure,  the  pavement  structure  (layer  thicknesses)  must 
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be  known  and  typically  moduli  can  be  backcalculated  for  only  two  or 
three  layers.  Therefore,  it  is  evident  that  a  fast,  economical,  and 
totally  nondestructive  method  for  determining  in  situ  modulus  profiles 
for  pavement  structures  is  needed. 

SASW  is  a  nondestructive  method  for  measuring  shear  wave  velocity 
and  determining  in  situ  elastic  moduli  of  soil  and  pavement  profiles  as 
well  as  the  thickness  of  each  layer.  SASW  is  based  on  the  generation 
and  detection  of  Rayleigh  waves  at  strains  below  0.001  percent  where 
the  elastic  properties  of  the  materials  are  independent  of  strain 
amplitude.  The  theory  of  elastic  waves  in  layered  media  is  utilized  to 
reduce  and  analyze  data  collected  in  the  field.  Practical  and 
theoretical  aspects  of  SASW  have  been  documented  (Nazarian  and  Stokoe, 
1985,  1986a).  SASW  has  been  under  continuous  development  at  the 
University  of  Texas  since  1980  and  at  the  Universities  of  Kentucky  and 
Michigan  since  1985.  A  review  of  the  historical  development  by 
Nazarian  (1984)  revealed  that  the  concepts  of  SASW  are  not  new, 
however,  civil  engineering  applications  were  limited  prior  to  the 
appearance  of  modern  spectral  analyzers  and  theoretically  based  methods 
of  determining  the  shear  wave  velocity  profile  from  collected  data 
(Nazarian,  1984).  The  use  of  SASW  for  nondestructive  pavement 
evaluation  is  gaining  in  popularity  and  a  number  of  case  studies  have 
been  documented  (Helsey,  Stokoe,  and  Meyer  1982;  Nazarian  and  Stokoe 
1983;  Nazarian,  Stokoe,  and  Hudson  1983;  Nazarian  and  Stokoe  1984; 
Drnevich,  et  al.  1985;  Nazarian  and  Stokoe  1986;  and  Hiltunen  1988). 
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B .  Purpose 

The  purpose  of  this  study  is  to  evaluate  the  applicability  of  SASW 
for  material  characterization  of  pavement  systems  in  situ  for  the 
determination  of  load  capacity  and  rehabilitation  requirements. 

C.  Scope 

The  efficiency  and  practical  feasibility  of  the  SASW  method  were 
evaluated  through  a  field  testing  program  designed  to  include  flexible, 
rigid,  and  composite  (flexible  over  rigid)  pavement  structures.  SASW 
field  testing  and  data  reduction  techniques  were  utilized  to  obtain 
Young's  modulus  pr.files  for  a  total  of  six  airfield  pavement  sites. 

For  comparison,  modulus  values  were  backcalculated  using  the  WES 
elastic  layer  method  and  FWD  surface  deflections.  SASW  moduli  were 
also  compared  to  laboratory  resilient  modulus  and  dynamic  modulus 
values  for  the  asphaltic  concrete  (AC)  and  Portland  cement  concrete 
(PCC)  materials,  respectively. 
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II.  SPECTRAL-ANALYSIS-OF-SURFACE-WAVES  (SASV)  METHOD 


A.  Wave  Propagation 

For  evaluation,  most  pavement  sites  can  be  approximated  by  a 
layered  half- space  since  the  depth  of  concern  is  usually  not  more  than 
20  ft.  A  pavement  system  can  be  modeled  by  a  number  of  elastic, 
isotropic,  homogeneous  layers.  Constant  properties  are  assumed  within 
each  layer.  If  waves  are  assumed  plane,  the  solution  of  wave  equations 
reduces  to  a  two-dimensional  problem. 

If  an  elastic  half-space  is  disturbed  by  a  vertical  Impact  on  the 
surface,  body  and  surface  waves  will  propagate  in  the  medium.  Body 
waves  propagate  radially  outward  in  the  medium  and  are  composed  of 
compression  waves  (P-waves)  and  shear  waves  (S-waves).  Particle  motion 
is  parallel  to  the  direction  of  propagation  for  P-waves  and 
perpendicular  to  the  direction  of  propagation  for  S-waves.  Surface 
waves  resulting  from  a  vertical  excitation  are  primarily  Rayleigh  waves 
(R- waves)  which  propagate  away  from  the  impact  along  a  cylindrical 
wavefront  near  the  surface  with  particle  motion  near  the  surface 
forming  a  retrograde  ellipse.  The  amplitude  of  the  wave  decays 
exponentially  with  depth,  such  that  at  a  depth  of  one  wavelength, 
vertical  and  horizontal  particle  motion  are  only  about  30  and  10 
percent  of  surface  values.  It  has  been  shown  (Miller  and  Pursey,  1955) 
that  approximately  67  percent  of  the  input  energy  propagates  in  the 
form  of  R-waves  while  P-  and  S-waves  carry  about  26  and  7  percent  of 
the  energy,  respectively.  As  the  waves  propagate,  the  Increasing 
volume  of  medium  encountered  causes  wave  energy  to  dissipate  or 
geometrical  damping.  While  body  waves  are  faster  than  surface  waves. 
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they  also  attenuate  much  faster  due  to  geometrical  damping.  Near  the 
surface,  body  waves  attenuate  proportional  to  1/r  (where  r  is  the 
distance  from  the  source,  while  R-waves  attenuate  proportional  to  l/Tr. 
Thus,  in  a  relatively  short  distance  from  the  source,  most  of  the 
energy  at  the  surface  is  R-wave  energy.  This,  in  addition  to  the  fact 
that  most  of  the  energy  from  an  impact  goes  into  R-wave  velocity,  is 
beneficial  to  the  SASW  method. 

The  velocity  of  propagation  of  the  different  waves  are  related  by 
Poisson's  ratio.  P-waves  propagate  faster  than  S-waves  and  the  ratio 
of  P-  to  S-wave  velocity  increases  from  1.4  to  infinity  as  the 
Poisson's  ratio  increases  from  0.0  to  0.5.  R-waves  are  the  slowest 
with  the  ratio  of  R-  to  S-waves  varying  from  0.86  to  0.95  for  Poisson's 
ratios  of  0.0  and  0.5,  respectively.  The  velocity  of  propagation  is  a 
direct  indicator  of  the  stiffness  of  a  material  and  once  the  shear  or 
compression  wave  velocity  is  known,  the  Young's  modulus  and  shear 
modulus  can  be  determined  as: 

G  - 

E  -  2G(1  +  n)  ~  2pVg2(i  +  (11-2) 

or 

E  -  pVp2[(l  +  ^)(l  .  2#i)/(l  -  n)]  (II-3) 

where , 

G  -  shear  modulus, 

E  -  Young's  modulus, 
p  -  mass  density, 

Vg  -  shear  wave  velocity. 
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H  -  Poisson's  ratio,  and 
Vp  -  compression  wave  velocity. 

R-wave  velocity  will  be  constant  and  independent  of  frequency  in  a 
homogeneous,  elastic  half -space.  For  a  case  such  as  a  pavement  system, 
the  stiffness  of  the  layers  varies  with  depth  and  the  R-wave  velocity 
will  vary  with  frequency.  Different  frequencies  will  propagate  at 
different  velocities  and  an  average  velocity  (phase  velocity)  for  the 
material  to  a  depth  of  1-  to  1-1/2  wavelengths  will  be  measured.  The 
variation  of  this  phase  velocity  with  wavelength  (frequency)  is  called 
dispersion. 

B.  Frequency  Domain  Analysis 

Analyses  in  the  frequency  domain  are  not  new,  however  until 
recently,  algorithms  such  as  the  Fourier  transform  required  a 
prohibitive  amount  of  computation  time  and  were  not  efficient  enough 
to  be  applied  in  the  field.  Development  of  microprocessors  and  the 
fast  Fourier  transform  (FFT)  have  greatly  extended  the  capability  to 
measure  and  analyze  dynamic  systems  in  the  frequency  domain.  The  FFT 
and  spectral  analyses  are  now  widely  used  in  field  and  laboratory 
testing  and  are  key  components  of  SASV  testing. 

SASV  measurements  are  made  on  two  channels  of  data.  A  number  of 
different  measurements  can  be  made  once  the  signals  have  been  Fourier 
transformed  (transformed  into  the  frequency  domain) .  These 
measurements,  spectral  analyses,  are  basically  statistical  operations 
performed  on  the  signals.  This  involves  the  correlation  of  the  input 
or  output  with  itself  or  the  correlation  of  the  input  with  the  output, 
assuming  that  the  object  under  test  is  a  linear  system. 
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Comparison  of  two  signals  is  accomplished  easily  in  the  frequency 
domain.  Mathematical  operations  can  be  performed  more  efficiently. 

For  example,  an  integration  in  the  time  domain  is  a  simple 
multiplication  in  the  frequency  domain.  An  important  characteristic  of 
most  frequency  domain  measurements  is  that  triggering  does  not  have  to 
be  synchronized.  This  makes  it  possible  to  obtain  enhanced  records  by 
averaging  several  signals.  For  repeatable  signals,  the  average  of 
several  records  will  theoretically  be  free  of  random  background  noise. 
The  spectral  analyses  functions  of  interest  in  SASU  testing  are 
described  below. 

Cross  power  spectrum.  Signals  from  two  channels  obtained 
simultaneously  can  be  compared  utilizing  the  cross  power  spectrum, 
Gvv(f)>  defined  as; 

G^y(f)  -  s/(f)  *  Sy(f)  (II-4) 

where , 

Sjj*(f)  -  the  complex  conjugate  of  the  linear  spectrum  of  the 
input , 

Sy(f)  -  linear  spectrum  of  the  output. 

The  linear  spectrum  is  the  Fourier  transform  of  the  signal  recorded  in 

the  time  domain.  G  (f)  is  a  complex  function  and  can  be  represented 

xy 

with  real  and  imaginary  components  or  with  a  magnitude  and  phase.  The 
magnitude  is  a  measure  of  the  mutual  power  between  the  two  signals. 

The  phase  is  the  relative  phase  between  the  signals  at  each  frequency 
in  the  measurement  bandwidth .  Because  the  phase  is  relative , 
synchronized  triggering  is  not  necessary.  The  cross  power  spectrum  is 
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the  key  measurement  in  SASW  testing  and  is  used  to  determine  the  travel 
times  of  surface  waves  at  different  frequencies. 

Auto  power  spectrum.  The  auto  power  spectrum,  is  defined  as; 

Gxx  -  Sx(f)  *  Sx*(f)  (II-5) 

where , 


S^Cf)  -  linear  spectrum, 

Sx  -  complex  conjugate  of  the  linear  spectrum 
The  magnitude  of  the  auto  power  spectrum  is  equal  to  the  magnitude  of 
the  linear  spectrum  squared  and  is  the  energy  at  each  frequency  in  the 
measurement  bandwidth.  Because  the  linear  spectrum  is  multiplied  by 
its  complex  conjugate,  the  auto  power  spectrum  is  a  real  valued 
function  and  does  not  contain  any  phase  information,  but  can  be  used  to 
identify  dominant  frequencies. 

Coherence .  The  coherence  function,  7  (f),  is  a  real  valued 
function  with  a  value  between  zero  and  one  corresponding  to  the  ratio 
of  the  output  power  caused  by  the  measured  input  to  the  total  measured 
output  defined  as: 


GxX<f)  * 


(II-6) 


where  G^,  G„„,  and  G„  are  previously  defined  cross  and  auto  power 

yx  XA  yy 

spectra  averaged  in  the  frequency  domain  over  multiple  input  records. 
The  coherence  function  should  be  used  on  averaged  signals  since  it  will 
always  have  a  value  of  unity  for  one  signal  pair.  A  coherence  value 
close  to  one  indicates  good  correlation  between  the  input  and  the 
output  signals.  A  low  coherence  may  indicate  bad  correlation,  however, 
nonlinearity  of  the  system  and  low  frequency  resolution  can  also  result 
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in  poor  coherence .  The  coherence  function  is  used  in  SASW  testing  to 
help  identify  frequencies  providing  accurate  cross  power  spectrum  phase 
values . 

C.  Overview 

SASW  is  a  nondestructive  seismic  method  for  determining  in  situ 
elastic  moduli  of  pavement  systems.  To  conduct  SASW,  two  receivers  are 
attached  to  the  surface  and  a  hammer  is  used  to  deliver  a  transient 
impact  to  the  pavement.  The  impact  generates  energy  over  a  range  of 
frequencies.  Wavelengths  depend  on  the  stiffness  of  the  material  and 
frequency  of  the  wave.  The  energy  of  a  surface  wave  decays  rapidly 
with  depth,  therefore,  different  frequencies  sample  different  depths. 

The  signal  generated  by  the  impact  is  recorded  using  a  waveform 
analyzer  and  the  time-domain  record  from  each  receiver  is  transformed 
into  the  frequency  domain.  This  process  is  repeated  several  times  and 
the  results  averaged  to  enhance  the  quality  of  the  signal  by 
eliminating  the  effects  of  random  background  noise. 

Spectral  analyses  are  then  performed  to  obtain  the  cross  power  and 
coherence  functions.  The  coherence  value  is  used  to  monitor  the 
quality  of  data  during  testing  and  to  identify  poor  data  during 
analysis.  The  cross  power  spectrum  provides  the  relative  phase  shift 
at  each  frequency  from  which  the  travel  time  between  the  two  receivers 
can  be  determined. 

Because  the  frequencies  that  can  be  generated  by  a  single  source 
are  limited,  testing  is  performed  at  several  spacings.  Close  spacings 
sample  material  near  the  surface  and  the  source  must  generate  high 
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frequencies.  Larger  spacings  sample  deeper  layers  and  the  source  must 
generate  more  energy  in  the  lower  frequency  ranges. 

The  first  step  in  data  reduction  involves  converting  the  field 
data  into  a  dispersion  curve.  The  dispersion  curve  is  simply  a  plot  of 
the  surface  wave  velocity  versus  wavelength  and  can  be  developed  from 
phase  information  of  the  cross  power  spectrum.  The  R-wave  velocities 
determined  in  this  manner  are  not  actual  velocities,  but  apparent 
(phase)  velocities  since  the  stiffness  of  layers  near  the  surface  will 
affect  the  measurement  of  velocities  of  underlying  layers. 

Finally,  the  shear  wave  velocity  profile  must  be  obtained  from  the 
dispersion  curve.  This  process  is  called  inversion  and  consists  of 
determining  the  depth  of  each  layer  and  the  actual  shear  wave  velocity 
of  each  layer  from  the  apparent  R-wave  velocity  versus  wavelength 
information.  Once  the  shear  wave  velocities  are  known,  the  shear  and 
Young's  modulus  can  be  computed. 
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III .  DATA  COLLECTION 

A.  Test  Sites 

Test  sites  for  this  study  were  selected  based  on  pavement 
structure  and  the  availability  of  adequate  data  to  verify  pavement 
conditions  at  the  time  of  SASW  testing.  The  six  airfield  pavement 
sites  listed  in  Table  1  represent  a  wide  range  of  pavement  types, 
thicknesses,  and  subgrade  conditions.  These  sites  were  selected  in 
conjunction  with  a  separate  study  that  was  being  conducted  by  WES  to 
evaluate  various  NDT  devices  and,  for  reference  purposes,  site  numbers 
in  Table  1  correspond  to  the  previously  reported  values  (Bentsen,  Bush, 
and  Harrison  1988).  During  that  study,  a  test  pit  was  excavated  within 
each  site  at  a  reference  point  location  for  structure  verification,  in- 
place  strength  measurements,  and  sampling  of  materials  for  laboratory 
testing.  Field  and  laboratory  test  results,  extracted  from  the 
previous  report,  for  each  of  the  six  sites  included  in  this  study  are 
summarized  in  Tables  2-6.  Table  2  is  a  siimmary  of  test  pit  data 
including  layer  thicknesses  and  in-place  density,  water  content,  and 
strength  measurements.  In-place  strengths  are  reported  in  terms  of  the 
California  bearing  ratio  (CBR)  for  flexible  pavement  layers  and  the 
modulus  of  subgrade  reaction,  k,  (from  a  plate  bearing  test)  for  rigid 
pavement  layers.  Base  and  subgrade  material  classifications  are 
presented  in  Table  3.  Laboratory  resilient  moduli  for  the  AC  layers 
determined  using  procedures  described  in  ASTM  D  4123-82  (ASTM  1988)  are 
shown  in  Table  4  for  temperatures  of  77-  and  104-degrees  fahrenhelt. 
Flexural  strengths  determined  using  ASTM  C-78-84  (ASTM  1986)  and 
Young's  modulus  values  from  ASTM  C- 215 -85  (ASTM  1988)  for  the  PCC 
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TABLE  1.  GENERAL  SITE  DESCRIPTION 


Pavement 

Thickness 


Site 

Airport/Airfield 

in. 

3 

NAS  Pensacola,  FL 

10.0 

4 

NAS  Pensacola,  FL 

5.5 

8 

Birmingham  ANG/ 

7.0 

Birmingham 
Municipal,  AL 

7.0 

10 

Birmingham  ANG/ 

2.0 

Birmingham 
Municipal,  AL 

7.0 

11 

Sheppard  AFB 

Wichita  Falls,  TX 

21.0 

12 

Sheppard  AFB 

Wichita  Falls,  TX 

7.0 

Type 

Physical  Location 

PCC 

North  end  of  taxiway 
parallel  to  runway  1-19 

AC 

North  end  of  taxiway 
parallel  to  runway  1-19 

AC 

PCC 

Southwest  end  of  ANG 
aircraft  apron 

AC 

PCC 

Taxiway  E  between  runway 
18-36  and  holding  apron 

PCC 

Alert  apron  at  North  end 
of  airbase 

AC 

Acces'-  -axiway  North  of 
cor.trol  tower  between 
taxiway  E  and  school 
apron 
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TABLE  2.  TEST  PIT  DATA  (THICKNESS,  DENSITY,  WATER  CONTENT,  CBR,  AND  K) 
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20.0  Base  143  147  3.8  .  3.1  7.0  117 

-  Subgrade  .  .  109  127  15.7  27.0  13 

.  96  ---  21.6  41.0  4.3 


TABLE  3.  SOIL  CLASSIFICATiai  TEST  RESUlis 
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TABLE  4.  LABORATORY  RESILIENT  MODULUS  OF  AC  SURFACE  LAYERS 


Sample 

Average  Instantaneous 

Thickness 

Resilient 

Modulus,  psi 

ite 

Lift 

in. 

77  Deg  F 

104  Deg  F 

4 

1 

2.50 

1,568,500 

619,759 

2 

2.54 

1,348,381 

349,607 

8 

1 

2.68 

903,572 

400,887 

2 

2.74 

596,460 

257,620 

10 

1 

2.02 

804,228 

356,960 

12 

1 

2.17 

669,266 

648,946 

2 

1.70 

1,230,766 

757,317 

3 

2.84 

1,177,419 

1,181,251 

TABLE  5.  LABORATORY  TEST  RESULTS,  PCC  LAYERS 


Site 

No. 

Dynamic 

Modulus 

psi 

Flexural 

Strength 

psi 

3 

5,556,000 

905 

8 

5,487,000 

820 

10 

6,794,000 

745 

11 

5,363,000 

510 
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TABLE  6.  CONDITION  SURVEY  RESULTS 


Site 

PCI 

Rating 

Prevalent  Distresses 

3 

82 

Very  Good 

Low-severity  patching  (large  and  small) 
Low- severity  joint  and  corner  spalls 

h 

97 

Excellent 

None 

8 

64 

Good 

Low- severity  block  cracking 

10 

68 

Good 

Low- severity  longitudinal  and  transverse 
cracking 

Medium-severity  joint  reflection  cracking 

11 

68 

Good 

High- severity  joint  seal  damage 

Low-  and  medium-severity  small  patches 

12 

76 

Very  Good 

Low- severity  longitudinal  and  transverse 
cracking 
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layers  are  presented  in  Table  5.  Results  of  a  condition  survey 
conducted  to  determine  the  pavement  condition  index  (PCI)  for  each  site 
are  summarized  in  Table  6.  The  PCI  is  a  numerical  indicator  based  on  a 
scale  of  0  to  100  and  is  determined  by  measuring  pavement  surface 
distress  that  reflects  the  surface  condition  of  the  pavement. 

B.  SASW  Testing  Equipment 

SASW  test  equipment  used  during  this  investigation  included 
recording  devices  for  monitoring,  recording,  and  performing  spectral 
analysis  on  waveforms,  sources  for  generating  waveforms,  receivers  for 
measuring  waveforms,  and  a  power  supply. 

A  Hewlett  Packard  (HP)  Model  3562A  Dynamic  Signal  Analyzer 
equipped  with  an  HP  Model  9153A  combined  10  megabyte  hard  disk  and  3- 
1/4  in.  floppy  disk  drive  was  used  in  this  study  for  collecting  and 
storing  SASW  data  records.  The  HP  3562A  is  a  versatile  2-channel 
device  with  a  sampling  rate  of  up  to  100  kHz  per  channel.  Spectral 
analyses  functions  such  as  the  FFT,  cross  power  spectrum,  auto  power 
spectrum,  and  coherence  can  be  performed  immediately,  thereby  allowing 
the  frequency  domain  functions  to  be  viewed  in  the  field. 

Impact  excitation  was  used  throughout  this  study.  The  Impact 
sources  were  an  array  of  hammers  including  a  4  oz.  ball  peen,  8  oz. 
ball  peen,  40  oz.  sledge,  and  8  lb.  sledge. 

Both  velocity  transducers  and  accelerometers  were  used  in  this 
study.  All  receivers  were  vertically  oriented.  Geosource  Model  PC-3 
velocity  transducers  were  used  for  low  frequency  measurements  (below 
1000  Hz).  PCB  Piezotronics ,  Inc.  Model  308B02  accelerometers  (resonant 
frequency:  32  kHz;  Sensitivity:  1000  mV/g) ,  powered  by  a  PCB  Model 
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483A  12  channel  power  supply,  were  used  for  high  frequency  measurements 
associated  with  pavement  surface  layers.  Coupling  between  the  pavement 
surface  and  receivers  was  accomplished  using  a  No.  2  plastilina 
sculpting  clay. 

The  power  supply  used  throughout  this  study  was  a  Honda  Model 
EM2200X  portable,  gasoline  generator  with  a  maximum  power  of  2.2  leva. 

An  Omega  Engineering,  Inc.,  Omegascope  infrared  temperature  gun  was 
used  for  determining  the  surface  temperature  of  asphalt  surface  layers 
at  the  time  of  testing. 

C.  SASW  Testing  Procedures 

SASW  tests  were  conducted  at  one  location  within  each  test  site 
corresponding  to  the  reference  point  location  referred  to  previously. 
The  common  receivers  midpoint  (CRMP)  geometry  test  procedure  shown 
schematically  in  Figure  1  was  used  for  all  SASW  testing  with  the 
centerline  located  at  the  reference  point.  Distances  of  0.25- ,  0.5- , 
1.0- ,  2.0- ,  4.0- ,  and  8.0-ft.  were  marked  off  on  either  side  of  the 
Imaginary  centerline  corresponding  to  spaclngs  of  0.5- ,  1.0- ,  2.0- , 

4.0- ,  8.0- ,  and  16.0-ft.  which  were  typically  used.  Data  were 
collec'  '  for  additional  spaclngs  only  when  deemed  necessary  based  on 
visual  oL  ervations  of  the  data  records  during  testing. 

It  should  be  noted  in  Figure  1  that  two  tests  are  Indicated  for 
each  spacing  with  the  source  on  opposite  sides  of  the  imaginary 
centerline.  These  are  termed  forward  and  reverse  profiles  and  can  be 
averaged  to  minimize  the  effect  of  any  internal  phase  shift  between 
receivers . 
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Figure  1.  Schematic  of  SASW  test  arrangement  (after 
Nazar ian,  1984) 
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The  purpose  of  multiple  spacings  is  to  insure  that  usable  data  are 
obtained  over  the  complete  range  of  frequencies  desired.  It  has  been 
shovm  (Nazarian  and  Stokoe  1983)  that  the  use  of  these  spacings  reduces 
scatter  in  the  data  due  to  overlap  between  consecutive  spacings.  Close 
receiver  spacings  are  associated  with  high  frequencies  and  sample 
shallow  depths  or  near  surface  layers  while  large  receiver  spacings  are 
associated  with  lower  frequencies  and  sample  deeper  materials. 

The  steps  followed  for  data  collection  at  each  site  were  as 
follows : 

1.  Initial  HP  3562A  Setup  for  SASW.  To  prepare  the  HP  3562A  for 
SASW  testing,  the  following  measurement  parameters  were  specified: 
linear  resolution,  frequency  response,  a  uniform  window,  and  the 
request  for  five  averages.  By  setting  the  OV  REJ  to  ON,  no  data  were 
accepted  if  either  channel  over-ranged.  Manual  preview  was  turned  ON 
to  help  establish  appropriate  ranges  and  bandwidths  for  each  spacing. 

It  was  turned  off  when  acceptable  data  were  obtained  for  a  receiver 
spacing.  The  screen  display  was  set  to  show  the  phase  of  the  cross 
power  spectrum  (in  ±  180  degree  format  or  modulo  2n)  and  the  coherence. 
Finally,  the  device  was  set  to  trigger  off  of  channel  1  input.  While 
other  parameters  including  the  frequency  span,  voltage  ranges  for 
channels  1  and  2,  the  trigger  level,  and  the  trigger  delay,  were 
adjusted  during  testing  depending  on  the  particular  spacing  and  source, 
this  completes  the  initial  setup  of  the  analyzer. 

2.  Receiver  Selection.  A  combination  of  accelerometers  and 
velocity  transducers  were  used  for  testing  the  six  pavement  sites. 
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Accelerometers  were  typically  used  for  spacings  of  4  ft.  and  below. 
Velocity  transducers  were  generally  used  for  spacings  of  4  ft.  and 
above.  The  range  of  frequencies  generated  for  each  spacing  was  used  to 
determine  when  to  switch  receiver  types.  Accelerometers  were  generally 
used  down  to  frequencies  of  500  to  1000  Hz,  The  standard  procedure  was 
to  conduct  two  series  of  tests  using  both  types  of  receivers  at  the 
spacing  where  a  change  in  receiver  types  was  made. 

3.  Selection  of  Sources.  The  initial  selection  of  a  source  was 
made  based  on  the  spacing.  Typically,  the  4  oz.  ball  peen  was  selected 
for  spacings  of  2  ft.  or  less,  the  8  oz.  ball  peen  or  40  oz.  sledge  for 
spacings  of  4-  to  8-ft. ,  and  the  8  lb.  sledge  for  the  16  ft.  spacing. 
The  actual  frequencies  required  for  each  spacing  are  dependent  on  the 
stiffness  of  the  materials,  thus,  the  final  hammer  selection  was  made 
upon  visual  examination  of  the  first  data  records. 

4.  Collection  of  Data.  With  a  receiver  pair  in  place,  the  source 
was  positioned  on  one  side  or  the  other  (forward  profile)  such  that  the 
source  to  near  receiver  distance  was  equal  to  the  spacing.  The 
receiver  nearest  the  source  was  always  connected  to  channel  1  on  the  HP 
3562A.  Several  sharp,  crisp  blows  (no  bounce)  of  the  hammer  were 
normally  required  to  adjust  the  remaining  HP  3562A  measurement 
parameters.  The  frequency  span  was  typically  set  such  that  more  than 
three -fourths  of  the  frequency  range  contained  high  quality  (high 
coherence)  data.  The  voltage  ranges  were  set  as  low  as  possible 
without  overloading  to  obtain  the  best  resolution  in  the  analog  to 
digital  conversion.  With  the  manual  preview  function,  time  signals 
could  be  viewed  after  the  first  impact  to  determine  if  further 
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adjustments  were  necessary.  After  the  frequency  span  and  voltage 
ranges  were  established,  the  trigger  level  was  set  to  trigger  about 
halfway  up  the  first  major  pulse  on  channel  1.  A  pre- trigger  delay  was 
then  set  for  channels  1  and  2  to  insure  that  the  initial  portion  of  the 
signal  was  obtained.  The  pre -trigger  delay  was  set  at  approximately  10 
percent  of  the  total  time  record.  These  parameters  generally  had  to  be 
set  for  each  receiver  spacing. 

After  all  the  measurement  parameters  were  set,  the  manual  preview 
was  turned  off  and  signals  from  five  impacts  collected.  The  phase  and 
coherence  records  were  saved  on  disk.  The  file  naming  convention  used 
was  a  three  digit  numer' cal  value  for  each  pair  of  records  preceded  by 
a  "P"  or  "C"  for  phise  and  coherence.  Then,  without  relocating  the 
receivers,  the  same  test  was  performed  with  the  source  on  the  opposite 
side  of  the  receivers  (reverse  profile)  with  the  receiver  cables 
switched  on  the  HP  3562A.  Pertinent  testing  Information  was  recorded 
on  specially  designed  SASW  data  forms. 

Steps  3  through  5  were  repeated  for  each  receiver  spacing.  SASW 
field  data  obtained  following  these  procedures  are  presented  in 
Appendices  A-F  (located  immediately  following  the  main  text)  for  Sites 
3,  4,  8,  10,  11,  and  12  respectively. 
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IV.  ANALYSIS  OF  DATA 

SASW  data  reduction  was  accomplished  in  two  phases.  First,  the 
field  data  for  each  site  was  converted  to  a  dispersion  curve.  The 
second  phase  then  consisted  of  determining  the  shear  wave  velocity 
profile  by  inverting  the  dispersion  curve. 

A.  Determination  of  Dispersion  Curve 

The  phase  of  the  cross  power  spectrum  and  coherence  records 
obtained  from  spectral  analyses  were  both  used  during  the  construction 
of  dispersion  curves.  While  the  actual  dispersion  curve  was  computed 
from  the  phase  information,  the  coherence  record  was  used  to  help 
identify  frequency  ranges  containing  "good"  data.  Nazarian  and  Stokoe 
(1986)  showed  that  frequencies  having  a  coherence  value  greater  than 
0.90  provide  useful  data.  While  this  criterion  was  not  rigidly 
followed,  the  coherence  was  used  extensively  to  evaluate  the  quality  of 
data.  The  phase  of  the  cross  power  spectrum,  with  "poor"  data  ranges 
removed  was  used  to  calculate  the  dispersion  curve. 

For  a  particular  frequency,  f,  a  phase  difference  of  360  degrees 
is  associated  with  a  travel  time  between  receivers  of  one  period  (1/f). 
Knowing  the  phase  for  a  particular  frequency  (^) ,  the  travel  time,  t, 
can  be  computed  as : 

t  - .  (IV-1) 

360  (f) 

The  phase  velocity  for  a  given  frequency  can  then  be  obtained  as 


X 


where  X  is  the  distance  between  receivers  and  the  corresponding 


23 


wavelength  is: 


Repeating  this  procedure  for  all  frequencies  yields  the  dispersion 
curve  for  the  given  receiver  spacing  and  source  location.  This  curve 
is  then  filtered  to  remove  all  wavelengths  greater  than  twice  the 
receiver  spacing.  These  wavelengths  are  too  long  for  a  given  receiver 
spacing  and  may  not  have  traveled  a  sufficient  distance  to  adequately 
s£imple  a  depth  proportional  to  the  wavelength.  The  procedure  is 
repeated  for  all  spacings  and  source  locations  and  the  data  combined  to 
construct  the  total  dispersion  curve  to  represent  the  site. 

The  computer  program  SASV,  containing  data  analysis  routines 
presented  by  Hiltunen  (1988),  was  used  to  reduce  the  field  data 
collected  during  this  study  and  construct  the  dispersion  curves.  The 
computer  program  SASW  requires  the  existence  of  a  pre-established  data 
file  identifying  data  records  to  be  included  in  the  dispersion  curve 
and  filtering  Information  for  the  removal  of  bad  data.  Data  files 
created  for  each  site,  containing  the  following  information  for  all 
spacings  and  source  locations,  are  included  in  Appendices  A-F: 

Line  1:  Name  of  data  file  containing  phase  information 
Line  2:  Cutoff  frequency  (the  maximum  frequency  for  which  usable 
data  are  available  for  a  particular  record,  determined  by 
visual  examination  of  phase  and  coherence  records) 

Line  3:  Number  of  poor  data  ranges  to  be  removed  before 

constructing  dispersion  curve  (determined  by  visual 


2k 


examination  of  phase  and  coherence  records) 

Lines  4^^  through  4j^  (N  -  number  of  poor  data  ranges);  Range 

containing  poor  data  that  are  to  be  removed  from  phase 
information  (for  each  range  of  poor  data,  the  starting 
frequency,  ending  frequency,  and  number  of  phase  wraps 
occurring  prior  to  the  ending  frequency  must  be 
specified) . 

Line  5:  Receiver  spacing,  ft. 

The  computer  program  SASW  uses  the  above  information  and  performs  the 
following  steps  to  construct  the  dispersion  curve  for  a  site: 

a.  Converts  raw  phase  data  from  ±  180  degree  format  to  a 
continuous  format  and  eliminates  frequencies  containing  poor 
quality  data, 

b.  Redigitizes  unwrapped  phase  vs.  frequency  data, 

c.  Generates  dispersion  curve  data  from  converted  phase  data, 

d.  Run  dispersion  cur'/e  data  through  wavelength/receiver  spacing 
filter  to  remove  all  wavelengths  greater  than  twice  the 
spacing, 

e.  Generates  an  average  dispersion  curve, 

f.  Compacts  the  dispersion  curve  data  (reduces  files  to  a  more 
manageable  size  by  averaging  consecutive  data  points  when  the 
frequencies  are  within  two  percent  of  each  other  and  the 
velocities  are  within  four  percent  of  each  other) . 

Final  dispersion  curves  are  shown  in  Figures  2-7  for  sites  3,  4, 
8,  10,  11,  and  12  respectively.  The  transition  between  the  high 
velocity  portion  of  the  dispersion  curve,  corresponding  to  the  surface 
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Figure  3.  Ejg^erlinental  di^iersion  csurve  for  site  4  (5.5"  AC/13.5”  base/silty 
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Figure  4.  Eb^jerimental  dispersion  curve  for  site  8  (6.5”  AQ/7”  PCx:/sandy 
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Figure  6.  E}^)eriinental  di^)erslon  cxirve  for  site  11  (21"  POC/6"  base/clayey 


PHASE  VELOCITY,  EPS 


31 


Figure  7.  Experimental  di^jersion  curve  for  site  12  (7”  AC/20''  base/sandy 


layer  and  the  lower  velocity  portion  of  the  curve,  corresponding  to  the 
subgrade  is  much  smoother  and  more  continuous  for  the  flexible 
pavements  (sites  4  and  12)  than  for  the  rigid  and  composite  pavements 
(sites  3,  8,  10,  and  11).  This  is  due,  in  part,  to  the  large  contrast 
in  velocity  between  the  PCC  layers  and  the  underlying  base  or  subgrade 
material.  Many  of  the  fluctuations  in  the  rigid  and  composite 
dispersion  curves  may  be  the  result  of  reflected  waves,  which  could  be 
a  combination  of  reflected  surface  waves  and  direct  and  reflected  body 
waves . 

B .  Inversion 

Stiff  layers  near  the  surface  of  pavements  cause  a  shift  in  the 
measured  velocities  of  underlying  layers.  Inversion  is  the  process  of 
obtaining  the  true  shear  wave  velocity  profile  from  the  phase 
velocities  associated  with  the  dispersion  curve.  Inversion  of  the  six 
pavement  sites  included  in  this  study  was  accomplished  using  the 
computer  program  INVERT  developed  by  Dr.  Soheil  Nazarian  of  the 
University  of  Texas  (1984).  INVERT  is  an  interactive  computer  program 
developed  to  compute  a  theoretical  dispersion  curve  from  an  assumed 
shear  wave  velocity  profile  using  a  modified  version  of  the  Thomson 
(1950) -Haskell  (1953)  matrix  formulation  for  layered  elastic  media,  as 
developed  by  Thrower  (1965).  Some  minor  modifications  have  been  made 
by  Dr.  Dennis  Hiltunen  of  Penn  State  University  (1988)  for 
implementation  on  a  personal  computer  (PC) . 

For  inversion,  each  site  was  divided  into  a  number  of  layers  based 
on  the  known  pavement  structure.  In  some  cases,  surface  and  base 
course  layers  were  subdivided  into  two  or  three  thinner  layers.  For 
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each  site,  a  semi -infinite  layer  was  located  at  a  depth  of 
approximately  one -third  of  the  maximum  wavelength  found  in  the 
dispersion  curve.  The  subgrade  was  arbitrarily  divided  into  several 
layers  of  various  thicknesses.  A  mass  density  was  assigned  to  base 
course  and  subgrade  materials  based  on  field  density  measurements.  A 
unit  weight  of  150  pcf  was  assumed  for  AC  and  PCC  materials.  Poisson's 
ratioii  of  0.35,  0.15,  0.35,  and  0.4  were  assumed  for  AC,  PCC,  base 
course,  and  subgrade  materials  respectively.  An  initial  estimate  of 
the  shear  wave  velocity  of  each  layer  was  made.  INVERT  was  then  used 
to  compute  the  theoretical  dispersion  curve  for  the  assumed  velocities. 
Experimental  and  theoretical  curves  were  compared  and  the  assumed 
velocities  changed,  using  an  iterative  procedure,  until  the  two  curves 
matched  within  a  reasonable  tolerance. 

INVERT  searches  for  points  on  the  theoretical  dispersion  curve 
based  on  a  user  specified  range  of  frequencies  and  velocities  and 
number  of  Increments  to  be  used  in  the  range.  The  dispersion  curve  can 
be  more  accurately  defined  using  a  large  number  of  Increments  (finer 
mesh),  however,  the  computational  time  increases  significantly  as  the 
number  of  increments  increases.  For  this  study,  the  inversions  were 
performed  on  a  PC  (AT,  12.5  MHz)  and  a  choice  of  21  increments  was 
considered  about  optimum.  Also,  instead  of  choosing  a  wide  range  of 
frequencies  and  velocities,  the  dispersion  curve  was  divided  into  three 
or  four  sections  and  the  matching  procedure  performed  on  each  of  the 
narrow  ranges  Individually. 

The  final  velocity  profiles  providing  a  reasonable  match  between 
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theoretical  and  experimental  curves  were  assumed  to  be  the  true  shear 
wave  velocity  with  depth.  A  comparison  of  the  measured  and  theoretical 
dispersion  curves  for  sites  3,  4,  8,  10,  11,  and  12  are  shown  in 
Figures  8-13  respectively.  The  Young's  modulus  was  then  computed  for 
each  layer  using  equation  11-2. 

The  final  velocity  profiles  and  computed  layer  moduli  are 
summarized  in  Tables  7-12  for  the  six  sites.  During  inversion  of  the 
rigid  pavements  (sites  3  and  11),  it  was  observed  that  the  theoretical 
dispersion  curve  was  not  very  sensitive  to  the  velocity  of  the  base 
course  layers  and  the  reported  moduli  may  be  in  error.  This  was  due  to 
the  large  contrast  in  velocity  between  the  PCC  layer  and  the  underlying 
base  course  and  the  fact  that  the  base  course  layer  was  relatively 
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Figure  8.  Oonnpariscjn  of  measured  and  theoretical  diversion  curves  (after 
im4rsiOT)  for  site  3  (10”  PCXy4”  base/silty  sand) 
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Figure  9.  Ocnpeurison  of  measured  emd  theoretical  dispersion  curves  (after 
inversion)  for  site  4  (5.5"  AC/13.5"  base/silty  sand) 
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Figure  10.  Oonparison  of  measured  and  theoretical  di^^ersion  curves  (after 
inversiOT)  for  site  8  (6.5”  AC/7”  POC/sandy  clay) 
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Ocnparlson  of  measured  and  theoretical  dispersion  curves  (after 
inversion)  for  site  10  (2”  AC/7''  POC/14"  base/sanc^  clay) 
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Figure  12.  Oonpariscxi  of  measured  and  theoretical  dispersicai  curves  (after 
inversion)  for  site  11  (21''  PCX:/6"  base/clayey  sand) 


o 


o  o  o  o 

— (  C\2  CO 


Id  ‘HlOMSiaAVAi 

40 


Figure  13.  Oonparison  of  measured  and  theoretical  di^)ersion  curves  (after 
inversion)  for  site  12  (7”  AC/20"  base/san^  clay) 


TABLE  7.  YOUNG'S  MODULUS  PROFILE  FOR  SITE  3  AFTER  INVERSION 


Laver 

Material 

Type 

Thickness 

ft 

Shear 

Wave 

Velocity 

fps 

Poisson's 

Ratio 

Mass 

Density 

slugs 

Young’ s 
Modulus 
psi 

1 

PCC 

0.42 

9,517 

0.15 

4.7 

6,799,286 

2 

PCC 

0.42 

9,517 

0.15 

4.7 

6,799,286 

3 

Ba'^e 

0.33 

804 

0.35 

4.0 

48,481 

4 

Subgrade 

2.00 

357 

0.40 

3.8 

9,417 

5 

Subgrade 

2.00 

474 

0.40 

3.8 

16,601 

6 

Subgrade 

3.00 

800 

0.40 

3.8 

47,289 

7 

Subgrade 

1.00 

1,050 

0.40 

3.8 

81,463 

TABLE  8.  YOUNG'S  MODULUS  PROFILE  FOR  SITE  4  AFTER  INVERSION 


Laver 

Material 

Type 

Thickness 

ft 

Shear 

Wave 

Velocity 

fps 

Poisson' s 
Ratio 

Mass 

Density 

slugs 

Young ' s 
Modulus 
psi 

1 

AC 

0.23 

4,250 

0.35 

4.6 

1,557,891 

2 

AC 

0.23 

3,250 

0.35 

4.6 

911,016 

3 

Base 

0.56 

1,026 

0.35 

4.2 

82,898 

4 

Base 

0.56 

941 

0.35 

4.2 

69,732 

5 

Subgrade 

2.00 

429 

0.40 

3.9 

13,956 

6 

Subgrade 

1.00 

645 

0.40 

3.9 

31,549 

7 

Subgrade 

1.00 

1,037 

0.40 

3.9 

81 , 549 
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TABLE  9.  YOUNG'S  MODULUS  PROFILE  FOR  SITE  8  AFTER  INVERSION 


Shear 

Wave 

Mass 

Young' s 

Material 

Thickness 

Velocity 

Poisson’s 

Density 

Modulus 

Laver  Type 

ft 

fps 

Ratio 

slugs 

psi 

TABLE  10.  YOUNG'S  MODULUS  PROFILE  FOR  SITE  10  AFTER  INVERSION 


Laver 

Material 

Type 

Thickness 

ft 

Shear 

Wave 

Velocity 

fps 

Poisson's 

Ratio 

Mass 

Density 

slugs 

Young ' s 

Modulus 

psi 

1 

AC 

0.17 

3,600 

0.35 

4.7 

1,142, 100 

1 

2 

PCC 

0.58 

9,300 

0.15 

4.7 

6,492,756 

3 

Base 

0.59 

988 

0.35 

4.1 

75,041 

4 

Base 

0.59 

988 

0.35 

4.1 

75,041 

5 

Subgrade 

1.50 

424 

0.40 

4.2 

14,682 

6 

Subgrade 

2.00 

697 

0.40 

4.2 

39,674  1 

7 

Subgrade 

2.00 

955 

0.40 

4.2 

74,482 

8 

Subgrade 

2.00 

1,082 

0.40 

4.2 

95,609 

9 

Subgrade 

1.00 

1,200 

0.40 

4.2 

117,600 

1 

AC 

0.27 

3,960 

0.35 

2 

AC 

0.27 

3,960 

0.35 

3 

PCC 

0.58 

6,840 

0.15 

4 

Subgrade 

1.00 

428 

0.40 

5 

Subgrade 

2.00 

720 

0.40 

6 

Subgrade 

2.00 

1,000 

0.40 

7 

Subgrade 

2.00 

1,100 

0.40 

8 

Subgrade 

1.00 

1,200 

0.40 

■ 

1,381,941 

■ 

1,381,941 

.7 

3,512,169 

.2 

14,960 

.2 

42,336 

.2 

81,667 

.2 

98,817 

.2 

117,600 

A2 


TABLE  11.  YOUNG'S  MODULUS  PROFILE  FOR  SITE  11  AFTER  INVERSION 


Laver 

Material 

Type 

Thickness 

ft 

Shear 

Wave 

Velocity 

fps 

Poisson's 

Ratio 

Mass 

Density 

slues 

Young ' s 
Modulus 
psi 

1 

PCC 

0.88 

8,727 

0.15 

4.7 

5,717,329 

2 

PCC 

0.88 

8,622 

0.15 

4.7 

5,580,579 

3 

Base 

0.50 

810 

0.35 

4.4 

54,128 

4 

Subgrade 

2.00 

404 

0.40 

3.9 

12,377 

5 

Subgrade 

3.00 

364 

0.40 

3.9 

10,048 

6 

Subgrade 

2.00 

408 

0.40 

3.9 

12,624 

7 

Subgrade 

1.00 

900 

0.40 

3.9 

61,425 

TABLE  12.  YOUNG'S  MODULUS  PROFILE  FOR  SITE  12  AFTER  INVERSION 


Laver 

Material 

Type 

Thickness 

ft 

Shear 

Wave 

Velocity 

fps 

Poisson's 

Ratio 

Mass 

Density 

slues 

Young ' s 
Modulus 
psi 

1 

AC 

0.20 

4,367 

0.35 

4.6 

1,644,847 

2 

AC 

0.20 

2,939 

0.35 

4.6 

745,003 

3 

AC 

0.20 

2,221 

0.35 

4.6 

425,458 

4 

Base 

0.84 

747 

0.35 

4.6 

48,128 

5 

BAse 

0.84 

651 

0.35 

4.6 

36,553 

6 

Subgrade 

3.00 

559 

0.40 

3.9 

23,696 

7 

Subgrade 

3.00 

501 

0.40 

3.9 

19,034 

8 

Subgrade 

1.00 

650 

0.40 

3.9 

32,040 

V.  ELASTIC  LAYER  METHODOLOGY 


A.  Procedure 

A  study  conducted  at  WES  (Bush  1980)  has  shown  that  the  deflection 
basin  produced  by  applying  a  load  to  the  surface  of  a  pavement  with  an 
NDT  device  can  be  used  to  derive  the  Young's  modulus  of  pavement 
layers .  A  computer  optimization  routine  was  developed  to  determine  a 
set  of  modulus  values  that  provide  the  best  fit  between  a  measured 
deflection  basin  and  a  computed  deflection  basin  when  given  an  initial 
estimate  of  the  elastic  modulus  values  and  a  limiting  range  of  moduli. 

Initially,  a  set  of  modulus  values  is  assumed  and  the  theoretical 
deflection  is  computed  at  offsets  corresponding  to  the  measured 
deflections.  Subsequently,  each  modulus  is  varied  individually,  and  a 
new  set  of  deflections  is  computed  for  each  variation.  A  simplified 
description  of  the  Iterative  process  for  adjusting  the  modulus  values 
and  matching  the  deflection  basins  is  shown  in  Figure  14.  This 
illustration  is  for  one  deflection  and  one  layer.  For  multiple 
deflections  and  layers,  the  solution  is  obtained  by  developing  a  set  of 
equations  that  define  the  slope  and  intercept  for  each  deflection  and 
each  unknown  layer  modulus  as  follows: 

Log  (Deflectionj )  -  Aj^  +  Sj^CLog  Ej^) 

where 

A  -  intercept, 

S  -  slope, 

j  -  1,  2,  . . .  ND  [ND  -  number  of  deflections],  and 

i  -  1,  2,  ...  NL  [NL  -  number  of  layers  with  unknown  moduli] 
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LOG  DEFLECTION 


LOG  MODULUS 


Figure  14. 


Simplified  description  of  how  deflection  basins  are  matched 
in  WESDEF  (one  deflection  and  one  layer) 
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B .  WESDEF 


The  modulus  backcalculation  computer  program  WESDEF  was  developed 
(Van  Cauwelaert,  Alexander,  White,  and  Barker  1988)  by  incorporating 
the  elastic  layer  program  WESLEA  (Van  Cauwelaert,  Lequeux,  Delaunnois 
1986a  and  Van  Cauwelaert,  Lequeux,  Delaunnois  1986b)  as  a  subroutine 
for  the  optimization  routine. 

WESDEF  incorporates  a  layer  of  infinite  thickness  having  a  modulus 
of  elasticity  of  1,000, COO  psi  and  Poisson's  ratio  of  0.5  below  the 
subgrade  layer.  This  stiff  layer  is  at  a  depth  of  20  ft.  unless  soil 
profiles  indicate  the  need  for  some  other  representation  (i.e.,  shallow 
rock) . 

WESDEF  is  capable  of  handling  both  multiple  loads  and  variable 
Interface  conditions.  For  a  given  layer  (n)  and  underlying  layer 
(n+1) ,  the  interface  value  can  range  from  1  for  complete  adhesion 
between  the  layers  to  1,000  for  almost  frictionless  slip  between  the 
layers.  Frictionless  slip  is  usually  assumed  at  the  bottom  of  a  PCC 
layer  and  full  adhesion  is  generally  assumed  for  most  other  pavement 
materials . 

Normally  three  iterations  of  the  program  produce  a  set  of  modulus 
values  that  yield  a  deflection  basin  that  is  within  an  average  of  3 
percent  of  each  of  the  measured  deflections.  Limiting  iteration 
criteria  require  that  the  absolute  sum  of  the  percent  differences 
between  computed  and  measured  deflections  or  the  predicted  change  in 
modulus  values  be  less  than  10  percent. 

C.  Limitations  of  Elastic  Layer  Deflection  Based  Procedure 

Limitations  of  the  elastic  layer  approach  include  the  assumptions 
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of  elastic,  homogenous  layers  extending  horizontally  to  infinity.  For 
the  case  of  badly  cracked  pavements,  the  validity  of  these  basic 
assumptions  becomes  questionable.  Other  major  assumptions  are  that  the 
NOT  loading  frequency  is  such  that  the  dynamic  load  can  be  replaced  by 
an  equal  static  load  and  that  measured  dynamic  deflections  are  equal  to 
theoretical  static  deflections.  Also,  the  number  of  layers  for  which 
moduli  can  be  calculated  is  limited  to  four  layers  and  must  be 
restricted  to  two  or  three  layers  in  most  cases  to  minimize  non¬ 
uniqueness  of  the  solution.  Thus,  for  a  typical  pavement  structure 
having  a  surface,  base,  and  subgrade,  the  subgrade  modulus  must  be 
assumed  constant  over  the  entire  depth  to  the  rigid  boundary  (usually 
located  at  20  ft.).  Layer  thicknesses  must  be  known  for  the  procedure 
to  yield  accurate  results.  Therefore,  if  good  construction  records  are 
not  available,  coring  is  required  to  obtain  the  necessary  structural 
information,  ff  bedrock  is  present  within  20  ft.  of  the  surface,  an 
accurate  determination  of  the  depth  is  necessary. 

D.  Backcalculation  of  Moduli  for  Test  Sites 

The  computer  program  WESDEF  was  utilized  with  surface  deflections 
measured  with  a  Dynatest  FWD  to  backcalculate  layer  moduli  for  each  of 
the  pavements  included  in  this  study.  The  FWD  surface  deflection  data 
in  Table  13,  extracted  from  the  previous  report  discussed  earlier 
(Bentsen,  Bush,  and  Harrison),  were  collected  within  each  test  site  on 
the  reference  point  location. 

The  Dynatest  Model  8003  FWD  is  a  trailer -mounted,  impulse  load 
device  that  can  produce  a  load  between  1,500  and  24,000  lb.  The 
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single -pulse  transient  load  is  generated  by  a  weight  dropping  on  rubber 
pads  which  transmit  the  force  to  the  pavement  through  an  11.8-in.- 
dlameter  dense  rubber  plate  cushioned  with  a  thin  rubber  pad. 
Deflections  are  determined  using  seven  geophones  that  are  typically 
placed  1-ft.  apart  from  the  center  of  the  load,  but  the  outer  six  can 
be  varied  from  12  to  95  in.  away  from  the  plate.  The  falling  weight 
system  is  controlled  by  a  Hewlett-Packard  Inter  Personal  Computer  (IPC) 
and  can  produce  up  to  eight  loadings  selected  from  any  combination  of 
four  adjustable  drop  heights  in  a  1-2  minute  period.  Load  and 
deflection  data  are  recorded  on  paper  for  each  loading  and  can  be 
automatically  saved  to  a  magnetic  disk. 

The  modulus  values  predicted  by  WESDEF  for  sites  3,  4,  8,  10,  11, 
and  12  are  shown  on  the  computer  outputs  in  Figures  15-20.  For  each 
site,  an  attempt  was  made  to  compute  modulus  values  for  ell  layers  and 
this  was  successfully  done  for  sites  4,  8,  and  12.  WESDEF  was  unable 
to  predict  realistic  values  for  the  base  course  layers  beneath  the 
rigid  pavements  (sites  3  and  11)  and  a  composite  modulus  was  obtained 
for  the  combined  base  and  subgrade  layers.  WESDEF  was  unable  to 
predict  realistic  modulus  values  for  all  four  layers  at  site  10  and  a 
modulus  value  of  230,000  psi  was  assigned  for  the  thin  AC  surface  layer 
based  on  the  mean  pavement  temperature  at  the  time  of  testing  and 
Figure  21. 
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*************************WESDEF  OUTPUT  SUMMARY************************ 


NUMBER  OF  ITERATIONS  PERFORMED:  2 


DEFLECTIONS  COMPUTED  FOR  FINAL  MODULUS  VALUES 
********************************************* 


SENSOR 

MEASURED 

COMPUTED 

OFFSET 

DEFLECTION 

DEFLECTION 

POSITION 

IN. 

MILS 

MILS 

DIFFERENCE 

%  DIFF. 

******** 

****** 

ic'kic'k’kic'krk’k’k 

********** 

it'k'JtiHe'k'k'kit'k 

******* 

1 

3.0 

8.6 

8.7 

-0.1 

-0.9 

2 

12.0 

7.9 

8.1 

-0.2 

-2.1 

3 

24.0 

7.3 

7.2 

0.1 

2.0 

4 

36.0 

6.2 

6.2 

0.0 

0.4 

5 

48.0 

5.4 

5.2 

0.2 

3.4 

6 

60.0 

4.3 

4.3 

0.0 

-0.9 

7 

72.0 

3.5 

3.6 

-0.1 

-1.6 

ABSOLUTE  SUM. 

0.7 

11.3 

ARITHMETIC  SUM: 

0.3 

AVERAGE- 

0.1 

1.6 

FINAL  MODULUS  VALUES 
******************** 


LAYER 

MODULUS 

POISSON'S 

THICK. 

INTERFACE 

NO. 

MATERIAL  TYPE 

PSI 

RATIO 

IN. 

VALUE 

***** 

********************** 

********* 

********* 

****** 

********* 

1 

PCC 

8242664. 

0.15 

10.00 

1000. 

2 

SUBGRADE 

19264. 

0.40 

230.00 

1. 

3 

RIGID  BOUNDARY 

1000000. 

0.50 

1.00 

1. 

4 

RIGID  BOUNDARY 

1000000. 

0.50 

1.00 

1. 

5 

RIGID  BOUNDARY 

1000000. 

0.50 

SEMI -INF 

ABSOLUTE  SUM  OF  %  DIFF.  NOT  WITHIN  TOLERANCE 
CHANGE  IN  MODULUS  VALUES  WITHIN  TOLERANCE 


Figure  15.  WESDEF  output  summary  for  site  3  (10"  PCC/A" 
base/silty  sand) 


50 


*************************WESDEF  OUTPUT  SUMMARY************************ 


NUMBER  OF  ITERATIONS  PERFORMED;  3 


DEFLECTIONS  COMPUTED  FOR  FINAL  MODULUS  VALUES 
********************************************* 


SENSOR 

MEASURED 

COMPUTED 

OFFSET 

DEFLECTION 

DEFLECTION 

POSITION 

IN. 

MILS 

MILS 

DIFFERENCE 

%  DIFF. 

★**★★*★★ 

****** 

********** 

********** 

********** 

******* 

1 

3.0 

38.9 

39.6 

-0.7 

-1.8 

2 

12.0 

22.2 

21.8 

0.4 

1.8 

3 

24.0 

11.6 

11.7 

-0.1 

-0.8 

4 

36.0 

7.0 

7.3 

-0.3 

-4.8 

5 

48.0 

4.8 

5.0 

-0.2 

-4.1 

6 

60.0 

3.6 

3.6 

0.0 

1.0 

7 

72.0 

2.8 

2.6 

0.2 

7.0 

ABSOLUTE  SUM: 

2.0 

21.4 

ARITHMETIC  SUM: 

-1.8 

AVERAGE: 

0.3 

3.1 

FINAL  MODULUS  VALUES 
******************** 


LAYER 

MODULUS 

POISSON' 

S  THICK. 

INTERFACE 

NO. 

MATERIAL  TYPE 

PSI 

RATIO 

IN. 

VALUE 

***** 

********************** 

********* 

******** 

•*  ******  ' 

********* 

1 

AC 

206491. 

0.35 

5.50 

1. 

2 

BASE  OR  SUBBASE 

35844. 

0.35 

13.50 

1. 

3 

SUBGRAOE 

20491. 

0.40 

221.00 

1. 

4 

RIGID  BOUNDARY 

1000000. 

0.50 

1.00 

1. 

5 

RIGID  BOUNDARY 

1000000. 

0.50 

SEMI -INF 

REACHED  MAX  NO  OF  ITERATIONS 

ABSOLUTE  SUM  OF  %  DIFF.  NOT  WITHIN  TOLERANCE 

CHANGE  IN  MODULUS  VALUES  WITHIN  TOLERANCE 


Figure  16.  WESDEF  output  summary  for  site  4  (5.5"  AC/13.5" 
base/silty  sand) 
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*************************WESDEF  OUTPUT  SUMMARY************************ 


NUMBER  OF  ITERATIONS  PERFORMED:  3 


DEFLECTIONS  COMPUTED  FOR  FINAL  MODULUS  VALUES 
********************************************* 


SENSOR  MEASURED  COMPUTED 
OFFSET  DEFLECTION  DEFLECTION 


POSITION 

IN. 

MILS 

MILS 

DIFFERENCE 

%  DIFF. 

******** 

****** 

**********  ********** 

********** 

******* 

1 

3.0 

21.0 

21.8 

-0.8 

-4.0 

2 

12.0 

10.6 

10.7 

-0.1 

-1.2 

3 

24.0 

8.1 

7.5 

0.6 

7.3 

4 

36.0 

5.4 

5.0 

0.4 

7.0 

5 

48.0 

3.3 

3.4 

-0.1 

-2.1 

6 

60.0 

2.2 

2.3 

-0.1 

-5.5 

7 

72.0 

1.5 

1.7 

-0.2 

-11.4 

ABSOLUTE  SUM 

2.3 

38.5 

ARITHMETIC  SUM 

-10.0 

AVERAGE 

0.3 

5.5 

FINAL  MODULUS  VALUES 
******************** 

LAYER 

MODULUS 

POISSON' 

S  THICK. 

INTERFACE 

NO. 

MATERIAL  TYPE 

PSI 

RATIO 

IN. 

VALUE 

***** 

********************** 

********* 

******** 

•*  ****** 

********* 

1 

AC 

114603. 

0.35 

6.50 

1. 

2 

PCC 

942662 . 

0.15 

7.00 

1000. 

3 

SUBGRADE 

34479. 

0.40 

226.50 

1. 

4 

RIGID  BOUNDARY 

1000000. 

0.50 

1.00 

1. 

5 

RIGID  BOUNDARY 

1000000. 

0.50 

SEMI -INF 

REACHED  MAX  NO  OF  ITERATIONS 

ABSOLUTE  SUM  OF  %  DIFF.  NOT  WITHIN  TOLERANCE 

CHANGE  IN  MODULUS  VALUES  NOT  WITHIN  TOLERANCE 


Figure  17.  WESDEF  output  summary  for  site  8  (6.5"  AC/7" 
PCC/sandy  clay) 
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*************************WESDEF  OUTPUT  SUMMARY************************ 


NUMBER  OF  ITERATIONS  PERFORMED;  3 


DEFLECTIONS  COMPUTED  FOR  FINAL  MODULUS  VALUES 

********************ie********-<t*************** 

SENSOR 

MEASURED 

COMPUTED 

OFFSET 

DEFLECTION 

DEFLECTION 

POSITION 

IN. 

MILS 

MILS 

DIFFERENCE 

%  DIFF 

******** 

****** 

**•*****•*-** 

********** 

********** 

******1 

1 

3.0 

11.7 

12.8 

-1.1 

-9.2 

2 

12.0 

10.3 

10.3 

0.0 

-0.1 

3 

24.0 

8.6 

8.4 

0.2 

2.7 

4 

36.0 

6.7 

6.5 

0.2 

2.9 

5 

48.0 

5,0 

4.9 

0.1 

1.9 

6 

60.0 

3.5 

3.6 

-0.1 

-2.9 

7 

72.0 

2.7 

2.7 

0.0 

-1.5 

ABSOLUTE  SUM; 

1.8 

21.2 

ARITHMETIC  SUM* 

-6.1 

AVERAGE 

0.3 

3.0 

FINAL  MODULUS  VALUES 
********************* 


LAYER 

MODULUS 

POISSON'S 

THICK. 

INTERFACE 

NO. 

MATERIAL  TYPE 

PSI 

RATIO 

IN. 

VALUE 

***** 

********************** 

********* 

********* 

****** 

********* 

1 

AC 

230000. 

0.35 

2.00 

1. 

2 

PCC 

6928886. 

0.15 

7.00 

1000. 

3 

BASE  OR  SUBBASE 

17183. 

0.35 

14.00 

1. 

4 

SUBGRADE 

27372. 

0.40 

217.00 

1. 

5 

RIGID  BOUNDARY 

1000000. 

0.50 

SEMI -INF 

REACHED  MAX  NO  OF  ITERATIONS 

ABSOLUTE  SUM  OF  %  DIFF.  NOT  WITHIN  TOLERANCE 

CHANGE  IN  MODULUS  VALUES  NOT  WITHIN  TOLERANCE 


Figure  18.  WESDEF  output  summary  for  site  10  (2"  AC/7 
PCC/14"  base/sandy  clay) 


*************************WESDEF  OUTPUT  SUMMARY************************ 


NUMBER  OF  ITERATIONS  PERFORMED:  1 


DEFLECTIONS  COMPUTED  FOR  FINAL  MODULUS  VALUES 
********************************************* 


SENSOR 

MEASURED 

COMPUTED 

OFFSET 

DEFLECTION 

DEFLECTION 

POSITION 

IN. 

MILS 

MILS 

DIFFERENCE 

%  DIFF. 

******** 

****** 

********** 

********** 

********** 

******* 

1 

3.0 

3.1 

3.2 

-0.1 

-3.8 

2 

12.0 

2.9 

3.0 

-0.1 

-4.6 

3 

24.0 

2.8 

2.9 

-0.1 

-3.8 

4 

36.0 

2.7 

2.8 

-0.1 

-2.4 

5 

48.0 

2.6 

2.6 

0.0 

-0.3 

6 

60.0 

2.2 

2.4 

-0.2 

-11.0 

7 

72.0 

2.2 

2.3 

-0.1 

-3.4 

ABSOLUTE  SUM: 

0.7 

29.3 

ARITHMETIC  SOM: 

-29.3 

AVERAGE: 

0.1 

4.2 

FINAL  MODULUS  VALUES 
******************** 

LAYER 

MODULUS 

POISSON'S 

THICK. 

INTERFACE 

NO. 

MATERIAL  TYPE 

PSI 

RATIO 

IN. 

VALUE 

***** 

********************** 

********* 

********* 

****** 

********* 

1 

PCC 

12498614. 

0.15 

21.00 

1000. 

2 

SUBGRADE 

15837. 

0.40 

219.00 

1. 

3 

RIGID  BOUNDARY 

1000000 . 

0.50 

1.00 

1. 

4 

RIGID  BOUNDARY 

1000000. 

0.50 

1.00 

1. 

5 

RIGID  BOUNDARY 

1000000 . 

0.50 

SEMI -INF 

ABSOLUTE  SUM  OF  %  DIFF.  NOT  WITHIN  TOLERANCE 
CHANGE  IN  MODULUS  VALUES  WITHIN  TOLERANCE 


Figure  19.  WESDEF  output  summary  for  site  11  (21"  PCC/6" 
base/clayey  sand) 
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*************************weSDEF  output  summary************************ 


NUMBER  OF  ITERATIONS  PERFORMED:  3 


DEFLECTIONS  COMPUTED  FOR  FINAL  MODULUS  VALUES 
********************************************* 


SENSOR 

MEASURED 

COMPUTED 

OFFSET 

DEFLECTION 

DEFLECTION 

POSITION 

IN. 

MILS 

MILS 

DIFFERENCE 

%  DIFF. 

******** 

****** 

********** 

********** 

********** 

******* 

1 

3.0 

49.6 

51.4 

-1.8 

-3.6 

2 

12.0 

29.4 

28.0 

1.4 

4.9 

3 

24.0 

14.9 

14.6 

0.3 

2.2 

4 

36.0 

8.2 

8.9 

-0.7 

-8.0 

5 

48.0 

5.5 

6.0 

-0.5 

-8.9 

6 

60.0 

4.2 

4.3 

-0.1 

-2.4 

7 

72.0 

3.8 

3.2 

0.6 

15.9 

ABSOLUTE  SUM: 

5.4 

45.8 

ARITHMETIC  SUM: 

0.0 

AVERAGE; 

0.8 

6.5 

FINAL  MODULUS  VALUES 
******************** 


LAYER 

MODULUS 

POISSON'S 

THICK, 

INTERFACE 

NO. 

MATERIAL  TYPE 

PSI 

RATIO 

IN. 

VALUE 

***** 

********************** 

********* 

********* 

****** 

********* 

1 

AC 

126030. 

0.35 

7.00 

1. 

2 

BASE  OR  SUBBASE 

20839. 

0.35 

20.00 

1. 

3 

SUBGRADE 

16872 . 

0.40 

213.00 

1. 

4 

RIGID  BOUNDARY 

1000000. 

0.50 

1.00 

1. 

5 

RIGID  BOUNDARY 

1000000. 

0.50 

SEMI -INF 

REACHED  MAX  NO  OF  ITERATIONS 

ABSOLUTE  SUM  OF  %  DIFF.  NOT  WITHIN  TOLERANCE 

CHANGE  IN  MODULUS  VALUES  WITHIN  TOLERANCE 


Figure  20.  WESDEF  output  summary  for  site  12  (7"  AC/20" 
base/sandy  clay) 
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VI.  SUMMARY  AMD  DISCUSSION 

A .  Sumaary 

Applicability  of  SASV  to  pavement  testing  has  been  demonstrated 
through  field  testing  and  material  characterization  of  six  airfield 
pavement  sites.  SASW  is  a  relatively  new  procedure  based  on  the 
generation  and  detection  of  surface  waves.  Using  the  theory  of  elastic 
waves  in  layered  media,  data  are  reduced  and  analyzed  (computer 
programs  SASU  and  INVERT)  to  determine  a  continuous  Young's  modulus 
profile  with  depth  for  a  pavement  structure  as  well  as  the  thickness  of 
each  layer.  SASW  testing  is  performed  at  strains  below  0.001  percent 
and  the  elastic  properties  of  the  materials  tested  can  be  assumed 
ii>..:.ependent  of  strain  amplitude.  The  average  time  required  for  field 
data  collection  was  51  minutes  per  test.  Reduction  and  analysis  of 
SASW  data  required  1  to  1-1/2  days  per  site  using  a  PC  (AT).  Analysis 
time  could  be  reduced  significantly  by  performing  the  inversions  on  a 
mainframe  computer  or  a  faster  PC  (386  machine) .  It  should  also  be 
noted  that  the  layering  for  inversion  was  based  on  actual  pavement 
thicknesses  and  no  attempt  was  being  made  to  determine  layer 
thicknesses.  If  the  thicknesses  would  not  have  been  available, 
computational  time  requirements  would  have  made  it  necessary  to  use  a 
faster  computer  to  obtain  results  similar  to  those  presented.  The 
shear  wave  velocity  profiles  presented  herein  could  be  further  refined 
by  subdividing  the  layers  and  continuing  the  inversion  process. 

Fluctuations  in  the  SASW  field  data  records  due  to  reflected  waves 
presented  some  difficulty  in  the  interpretation  of  data  for  the 
construction  of  dispersion  curves  for  the  sites  with  PCC  layers.  The 
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effects  of  these  fluctuations  could  be  seen  in  the  dispersion  curve, 
making  it  difficult  to  perform  the  inversion  process  which  requires 
visual  matching  of  the  experimental  and  theoretical  data.  The  effect 
of  reflected  waves  on  SASW  testing  has  been  addressed  by  Sheu,  Stokoe, 
and  Roesset  (1988).  Another  problem  encountered  with  rigid  pavements 
was  the  relative  insensitivity  of  the  inversion  process  to  the  velocity 
of  a  thin  base  course  layer  beneath  a  PCC  layer. 

For  comparison  with  SASW,  the  WES  elastic  layer  NDT  methodology 
has  also  been  presented.  The  WES  procedure  utilizes  surface  deflection 
measurements  from  a  loading  device,  such  as  the  FWD,  and  elastic  layer 
theory  to  determine  a  set  of  modulus  values  that  provides  the  best  fit 
between  the  measured  and  computed  deflections.  Data  analyses  are 
accomplished  using  an  optimization  routine  (WESDEF)  that  performs  the 
iterative  basin  matching  process.  With  this  procedure,  modulus  values 
could  not  be  computed  for  thin  base  course  layers  located  directly 
beneath  PCC  layers. 

SASW  and  WESDEF  Young's  modulus  profiles  for  sites  3,  4,  8,  10, 

11,  and  12  are  presented  graphically  in  Figures  22-27  and  in  tabular 
form  in  Table  14.  Laboratory  derived  moduli  for  the  AC  and  PCC  layers 
are  included  in  Table  14. 

B.  Comparison  of  Moduli  for  Different  Material  Types 

Soil  (subgrade.  base,  and  subbase).  Hardin  and  Dmevich  (1972) 
have  suggested  that  the  main  factors  affecting  moduli  measured  in  the 
laboratory  are  state  of  stress,  void  ratio,  and  strain  amplitude.  For 
in  situ  tests  using  SASW,  strain  amplitude  will  have  essentially  no 
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Figure  23.  SASW  and  WESDEF  Young's  nicxSulus  profiles  for  site  4  (5.5"  AC/13.5" 
base/silty  sand) 
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(Continued) 
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effect  on  the  in  situ  tests  because  the  measurements  are  performed  at 
very  low  strains.  Numerous  studies  have  shown  that  shear  moduli  are 
nearly  constant  up  to  a  shearing  strain  of  0.01  percent  with  only  a 
slight  decrease  in  the  strain  range  from  0.001  to  0.01  percent.  The 
modulus  determined  in  this  constant  range  is  referred  to  as  the  low- 
amplitude,  initial  tangent,  or  maximum  modulus.  As  the  strain 
increases  above  0.01  percent,  moduli  decrease  significantly  and  are 
normally  determined  in  the  laboratory  or  are  estimated  from  empirical 
correlations . 

To  illustrate  the  variation  in  Young's  modulus  with  normal  strain, 
Nazarian  and  Stokoe  (1986)  presented  the  resonant  column  data  in 
Figures  28  and  29  for  a  stiff  clay.  A  constant  modulus  is  observed 
below  strain  levels  of  0.001  in  Figure  28  and  the  low-amplitude  modulus 
increases  with  increasing  confining  pressure.  In  Figure  29,  the 
noirmalized  modulus  is  also  constant  below  0.001  percent  and  is  equal  to 
the  maximum  modulus.  When  normalized,  all  of  the  modulus -strain  curves 
are  nearly  independent  of  confining  pressure  and  if  a  such  a  curve  is 
available  for  a  material,  the  moduli  at  higher  strains  can  be 
determined  once  the  maximum  modulus  has  been  measured  (i.e.,  small 
strain  in  situ  measurements  such  as  SASW) . 

For  the  type  and  magnitude  of  loading  applied  by  the  FWD,  the 
strain  levels  are  not  known  and  the  modulus  values  from  WESDEF  may  not 
be  the  maximum  values.  Therefore,  direct  comparisons  between  SASW  and 
WESDEF  values  for  the  soil  materials  will  be  questionable  without 
laboratory  test  results  to  verify  behavior.  Of  more  importance  here 
will  be  the  relative  stiffnesses  indicated  by  each  procedure. 
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Normalized  Young's  Modulus 


Single  Amplitude  Axial  Strain,  < ,  percent 

Figure  28.  Variation  of  Young's  modulus  with  strain  amplitude  at 
different  confining  pressures  of  an  unsaturated  clay 
subgrade  (from  Nazar ian  and  Stokoe,  1986a) 


Single  Amplitude  Axial  Strain,  c  ,  percent 


Figure  29.  Variation  of  normalized  Young's  modulus  with  strain 

amplitude  of  an  unsaturated  clay  subgrade  (from  Nazarian 
and  Stokoe,  1986b) 
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Base  course  moduli  from  SASW  were  typically  higher  than  WESDEF 
values  by  a  factor  of  at  least  two.  In  general,  SASW  subgrade  modulus 
values  were  lower  than  WESDEF  values  near  the  top  of  the  subgrade . 
Variation  of  SASW  moduli  with  depth  in  the  subgrade  was  different  for 
each  of  the  three  subgrades  evaluated.  The  sandy  subgrade  at  Pensacola 
(sites  3  and  4)  exhibited  a  gradual  increase  in  stiffness  from 
approximately  10,000  psi  to  values  of  50,000-  to  80,000-psi  at  depths 
of  6-  to  8-ft.  This  type  of  increase  in  stiffness  with  increasing 
confining  pressure  is  a  reasonable  behavior  for  a  sand.  WESDEF 
predicted  a  subgrade  modulus  of  approximately  20,000  psi  to  represent 
the  material  down  to  a  depth  of  20  ft.  This  is  approximately  one 
fourth  of  the  SASW  value  observed  at  a  depth  of  8  ft.  and  the 
difference  may  be  due  to  the  strain  amplitude  of  measurement.  Close 
agreement  between  the  subgrade  modulus  values  obtained  for  sites  3  and 
4  (located  within  2000  ft.  of  each  other)  was  observed  for  both 
procedures.  The  subgrade  at  Birmingham  (sites  8  and  10)  increased 
sharply  from  about  15,000  psi  to  values  of  100,000  psi  at  6-  to  8-ft. 
depths,  possibly  indicating  the  presence  of  bedrock  somewhere  near  the 
surface.  The  WESDEF  subgrade  modulus  values  were  also  higher  at  the 
two  Birmingham  sites  (approximately  30,000  psi).  Again,  good  agreement 
was  obtained  between  the  subgrade  moduli  computed  for  the  two  sites 
which  were  located  across  the  airfield  from  each  other.  SASW  subgrade 
modulus  values  were  relatively  constant  to  depths  of  8-  to  10-ft.  in 
the  sandy  clay  material  at  Wichita  Falls  (sites  11  and  12).  WESDEF 
modulus  values  were  approximately  the  same  for  the  two  sites,  located 
across  the  airfield  from  each  other.  SASW  subgrade  modulus  values  were 
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much  lower  for  site  11  (12,000  psi)  than  for  site  12  (20,000  psi). 

Asphaltic  concrete.  Variables  influencing  the  elastic  modulus  of 
asphaltic  concrete  include  temperature,  mix  properties,  age,  strain 
amplitude,  number  of  load  applications,  and  frequency  of  loading.  The 
two  main  factors  are  temperature  and  loading  frequency.  The  modulus 
generally  increases  with  decreasing  temperature  and  increasing 
frequency.  With  time,  an  aging  affect  tends  to  stiffen  the  asphaltic 
material  and  increase  the  modulus.  The  modulus  will,  however,  decrease 
with  increasing  strain  amplitude  and  load  repetitions  or  fatigue. 

Strain  amplitude  of  measurements  appeared  to  have  a  significant 
influence  on  the  modulus  of  the  AC  materials.  The  low  amplitude  SASW 
modulus  values  were  consistently  higher  than  the  WESDEF  values  (ranging 
from  three  to  ten  times  the  WESDEF  values).  The  SASW  values  seemed  to 
agree  more  closely  with  the  resilient  modulus  values  determined  at  77 
degrees  fahrenheit,  even  though  the  temperatures  at  the  time  of  testing 
were  closer  to  104  degrees.  SASW  modulus  values  for  sites  4  and  12 
decreased  significantly  with  depth  in  the  AC  layer.  Laboratory  values, 
which  are  highly  variable,  for  site  12  actually  showed  an  increase  in 
resilient  modulus  with  depth.  The  decrease  in  the  SASW  values  could 
not  be  explained,  but  may  be  due  to  the  influence  of  the  weaker  layer 
beneath  the  AC. 

Portland  cement  concrete .  In  general,  the  elastic  modulus  of 
Portland  cement  concrete  increases  with  increased  age,  rising  rapidly 
during  the  first  few  months  and  continuing  for  up  to  three  years.  Two 
major  factors  affecting  the  modulus  are  the  water-cement  ratio  and  the 
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type  of  aggregate.  High  cement  factor  and  low  water-cement  ratio  lead 
to  higher  elastic  moduli.  A  higher  modulus  can  be  expected  for 
concrete  made  with  stiffer  aggregate.  Many  mixing  and  curing  variables 
affect  the  strength  of  concrete  and  generally  will  influence  the 
modulus  in  a  similar  manner,  but  to  a  lesser  extent. 

PCC  moduli  from  SASW  appeared  more  reasonable  than  WESDEF  values 
in  several  cases.  WESDEF  predicted  values  over  8,000,000  psi  for  sites 
3  and  11,  which  would  normally  be  considered  unrealistically  high. 

SASW  modulus  values  agreed  better  with  laboratory  values  than  did 
WESDEF  values .  SASW  values  were  within  25  percent  of  laboratory  values 
for  three  out  of  four  sites  and  correlation  of  the  two  appears 
promising. 
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VII .  CONCLUSIONS 


The  major  conclusions  resulting  from  this  research  are  as  follows; 

1)  SASW  is  a  viable  method  for  in  situ  material  characterization 
of  pavement  systems.  In  situ  Young's  modulus  profiles  can  be 
obtained  for  flexible,  rigid,  and  composite  pavement 
structures.  The  ability  to  identify  layering  within  a 
pavement  structure  and  the  change  in  stiffness  within  layers 
has  many  applications  in  pavement  design  and  evaluation 
including  the  development  of  better  models  for  predicting  soil 
behavior  for  use  with  elastic  layer  or  more  advanced  finite 
element  routines . 

2)  SASU  data  reduction  can  be  accomplished  on  a  personal  computer 
(PC)  to  the  extent  of  identifying  relative  material 
stiffnesses  which  would  be  an  Immediate  benefit  to  the  current 
elastic  layer  deflection-based  procedure.  Determination  of 
pavement  layer  thicknesses  or  a  more  refined  Inversion  process 
will  require  a  mainframe  computer  or  a  fast  PC  (386  machine) . 

3)  PCC  modulus  values  from  SASW  appeared  reasonable  for  the  four 
rigid  and  composite  pavements  tested,  while  WESDEF  predicted 
unrealistically  high  values  for  two  of  the  sites.  SASW  values 
were  typically  within  25  percent  of  laboratory  dynamic  moduli 
and  the  existence  of  a  correlation  between  the  two  appears 
promising. 

A)  AC  moduli  determined  from  low  strain  amplitude  SASW 
measurements  were  three  to  ten  times  higher  than  the  values 
predicted  by  WESDEF.  Laboratory  resilient  modulus  values 
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determined  at  77  degrees  fahrenheit  agreed  more  closely  with  SASV 
moduli  than  did  104  degree  values,  which  should  be  closer  to  the 
actual  pavement  temperatures  measured  when  the  SASU  tests  were 
conducted. 

5)  The  effect  of  low  strain  amplitude  was  evident  in  the 
magnitude  of  SASW  base  and  subgrade  moduli.  SASW  moduli  were, 
in  some  cases,  typically  two  to  four  times  higher  than  WESDEF 
base  course  and  subgrade  values. 

6)  Low  amplitude  (maximum)  modulus  values  may  have  application 
ti--.  pavement  design  based  on  laboratory  test  results.  The 
ability  to  translate  between  laboratory  and  field  conditions 
using  SASU  appears  promising. 

7)  SASU  testing  and  data  reduction  techniques  are  somewhat  time 
consuming,  but  could  be  automated  to  greatly  increase  the 
efficiency. 
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VIII .  RECOMMENDATIONS 


This  study  has  shown  that  SASU  can  provide  certain  enhancements  to 
the  current  deflection-based  elastic  layer  procedure  and  may  have 
Important  applications  for  laboratory-based  pavement  design.  However, 
several  areas  requiring  further  research  have  been  Identified. 

Based  on  the  results  of  this  study,  the  following  recommendations  are 
offered: 

1)  Evaluate  methods  for  automating  data  collection  and  reduction 
techniques  to  Improve  the  efficiency  of  SASU  for  routine  pavement 
evaluation. 

2)  Evaluate  low  strain  effects  on  SASU  soil  moduli  by  comparing 
measured  values  to  laboratory  resilient  moduli. 

3)  Evaluate  methods  for  minimizing  the  effects  of  reflected  waves  on 
SASU  test  results. 

A)  Conduct  SASU  field  tests  on  pavements  with  bedrock  at  a  known 
depth  to  evaluate  the  ability  of  SASU  to  determine  its  presence 
and  approximate  depth. 

5)  Pursue  the  development  of  a  correlation  between  SASU  PCC  moduli 
and  the  laboratory  dynamic  modulus . 

6)  Further  Investigate  the  effects  of  low  strain  measurements  on  AC 
materials . 
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Figure  A2.  Phase  and  coherence  records  for  0.5  ft.  receiver 
spacing  at  Site  3 
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spacing  at  Site  3 
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Figure  A9.  Data  file  containing  the  names  of  phase  records,  cutoff 

frequencies ,  poor  data  ranges ,  and  receiver  spacings  used 
by  the  computer  program  SASW  in  constructing  the  dispersion 
curve  for  Site  3 
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Figure  Bl.  SASW  field  data  form  for  Site  A 
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Figure  B3.  Phase  and  coherence  records  for  1.0  ft.  receiver 
spacing  at  Site  4 
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Figure  B5.  Phase  and  coherence  records  for  4.0  ft.  receiver 
spacing  at  Site  4  (accelerometer  data) 


Figure  B6.  Phase  and  coherence  records  for  4.0  ft.  receiver 
spacing  at  Site  4  (velocity  transducer  data) 
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Figure  B9 .  Data  file  containing  Che  names  of  phase  records ,  cutoff 

frequencies,  poor  data  ranges,  and  receiver  spacings  used 
by  the  computer  program  SASW  in  constructing  the  dispersion 
curve  for  Site  4 
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APPENDIX  C:  SASW  TEST  DATA  AND  RESULTS  FOR  SITE  8 
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Figure  Cl.  SASU  field  data  form  for  Site  8 

C2 


Figure  C2.  Phase  and  coherence  records  for  0.5  ft.  receiver 
spacing  at  Sice  8 
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Figure  C3.  Phase  and  coherence  records  for  1.0  ft.  receiver 
spacing  at  Site  8 
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Figure  C4.  Phase  and  coherence  records  for  2.0  ft.  receiver 
spacing  at  Site  8 
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Figure  C6.  Phase  and  coherence  records  for  4.0  ft.  receiver 
spacing  at  Site  8  (velocity  transducer  data) 
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Figure  C8.  Phase  and  coherence  records  for  16.0  ft. 

receiver  spacing  at  Site  8  (250  Hz  bandwidth) 
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Figure  C9 .  Phase  and  coherence  records  for  16.0  ft. 

receiver  spacing  at  Site  8  (100  Hz  bandwidth) 
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Figure  CIO.  Data  file  containing  the  names  of  phase  records,  cutoff 

frequencies,  poor  data  ranges,  and  receiver  spacings  used 
by  the  computer  program  SASU  in  constructing  the 
dispersion  curve  for  Site  8 
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Figure  Dl.  SASW  field  data  form  for  Site  10 
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Figure  D2.  Phase  and  coherence  records  for  0.5  ft.  receiver 
spacing  at  Site  10 
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Figure  D3.  Phase  and  coherence  records  for  0.25  ft. 
receiver  spacing  at  Site  10 
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Figure  D8.  Phase  and  coherence  records  for  16.0  ft 
receiver  spacing  at  Site  10 
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Figure  DIO.  Data  file  containing  the  names  of  phase  records,  cutoff 

frequencies,  poor  data  ranges,  and  receiver  spacings  used 
by  the  computer  program  SASW  in  constructing  the 
dispersion  curve  for  Site  10 
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APPENDIX  E 


SASW  TEST  DATA  AND  RESULTS  FOR  SITE  11 
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Figure  El.  SASW  field  data  form  for  Site  11 
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Figure  E2.  Phase  and  coherence  records  for  0.5  ft.  receiver 
spacing  at  Site  11 
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Figure  E4.  Phase  and  coherence  records  for  2.0  ft.  receiver 
spacing  at  Site  11 


E5 


Unif 


E6 


IIMiSI 


lll'iU 


Fxd  Y  O  Hz 


C432 


CROSS  SPEC 


OSOvlD  Unlf 


mia 

I  ilM 
A  iMKara 

llllHIffiillS 


Fxd  Y  O 


Figure  E8.  Phase  and  coherence  records  for  16.0  ft. 

receiver  spacing  at  Site  11  (500  Hz  bandwidth) 
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receiver  spacing  at  Site  11  (250  Hz  bandwidth) 
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Figure  ElO.  Data  file  containing  the  names  of  phase  records,  cutoff 

frequencies,  poor  data  ranges,  and  receiver  spacings  used 
by  the  computer  program  SASW  in  constructing  the 
dispersion  curve  for  Site  11 
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APPENDIX  F 


SASW  TEST  DATA  AND  RESULTS  FOR  SITE  12 
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Figure  FI.  SASW  field  data  form  for  Site  12 
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Figure  F5.  Phase  and  coherence  records  for  4.0  ft.  receiver 
spacing  at  Site  12 
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Figure  F8.  Data  file  containing  the  names  of  phase  records,  cutoff 

frequencies,  poor  data  ranges,  and  receiver  spacings  used 
by  the  computer  program  SASW  In  constructing  the  dispersion 
curve  for  Site  12 
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